rm(list=ls())

source("bacom.R")


alpha1.hat <- numeric(0)

for (alpha in 0.1*(3:8)) {

  data.file.name <- paste("data/simData/simData_alpha_", alpha, ".RData", sep="")
  load(data.file.name)
  a <- bacom(simData$chrID, simData$copyNumber, simData$alleleA, simData$alleleB,
                 simData$isGenotypeAB, simData$segMean)  
  alpha1.hat <- c(alpha1.hat, a)
  
}



alpha2.hat <- numeric(0)

for (alpha in 0.1*(3:8)) {

  data.file.name <- paste("data/simData/simData_hemi_only_alpha_", alpha, ".RData", sep="")
  load(data.file.name)
  a <- bacom(simData$chrID, simData$copyNumber, simData$alleleA, simData$alleleB,
                 simData$isGenotypeAB, simData$segMean)  
  alpha2.hat <- c(alpha2.hat, a)
  
}


cat("alpha", "\t Dataset2", "\t Dataset1", "\n")
for (i in (1:6)) {

  cat((i+2)*0.1, "\t", round(alpha1.hat[i], 3), "\t", round(alpha2.hat[i], 3), "\n")

}
